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Enterovirus B81 (EV-B81) is a newly identified serotype within the species enterovirus B (EV-B). To date, 
only eight nucleotide sequences of EV-B81 have been published and only one ftiU-length genome sequence 
(the prototype strain) has been made available in the GenBank database. Here, we report the full-length 
genome sequences of two EV-B81 strains isolated in the Tibet Autonomous Region of China during acute 
flaccid paralysis surveillance activities, and we also conducted an antibody seroprevalence study in two 
prefectures of Tibet. The sequence comparison and phylogenetic dendrogram analysis revealed high 
variability among the global EV-B81 strains and frequent intertypic recombination in the non-structural 
protein region of EV-B serotypes, suggesting high genetic diversity of EV-B81. However, low positive rates 
and low titers of neutralizing antibodies against E V-B8 1 were detected. Nearly 68% of children under the age 
of five had no neutralizing antibodies against EV-B8 1 . Hence, the extent of transmission and the exposure of 
the population to this EV type are very limited. Although little is known about the biological and pathogenic 
properties of E V-B8 1 because of few research in this field owing to the limited number of isolates, our study 
provides basic information for further studies of EV-B81. 



Enteroviruses (EVs) belong to the family Picornaviridae within the new order Picornavirales. The genome of 
the small non-enveloped viruses is a single-stranded, positive-sense RNA molecule of approximately 7500 
nucleotides consisting of a single open reading frame flanked by 5' and 3' untranslated regions (UTRs). A 
single polyprotein translated from the RNA strand is first cleaved into three polyprotein precursors: PI, P2, and 
P3. PI is processed to yield four structural proteins: Viral protein 1-4 (VP1-VP4); P2 and P3 are precursors of the 
nonstructural proteins 2A-2C and 3A-3D, respectively. 

Most EV infections are asymptomatic or cause only mUd symptoms. However, some EVs can also cause a 
broad spectrum of other clinical illnesses, including acute flaccid paralysis (AFP); acute hemorrhagic conjunc- 
tivitis; encephalitis; aseptic meningitis; and hand, foot, and mouth disease' ''. Human EVs now comprise more 
than 110 serotypes, which are currently classified into four species, EV-A to EV-D, on the basis of their molecular 
and biological characteristics^. Although generally reliable, the neutralization test-a traditional method for EV 
typing-has gradually been replaced because it is labor-intensive and time-consuming; furthermore, it may fail to 
identify an isolate because of aggregation of virus particles, antigenic drift, or the presence of multiple viruses in 
the specimen''. Current EV classification is based on the high nucleotide sequence divergence within the VPl 
capsid-coding region, which has been shown to correspond with serotype neutralization*''. According to the 
recommended specific criteria for the interpretation of VPl sequence data, EVs are classified into the same type if 
they have more than 75% nucleotide similarity (85% amino acid similarity) and into different types if they have 
less than 70% nucleotide similarity in this region. While 70-75% nucleotide similarity in VPl region has been 
considered as a "grey zone" of molecular typing of EVs, and in this instance, additional information such as 
complete PI region sequences nucleotide similarity or neutralization profile should be obtained to decide the EV 
serotype'. A large number of new EV types have been discovered after molecular typing methods became 
available. 

Enterovirus B81 (EV-B81) is a newly identified serotype within the species EV-B. The prototype strain (USA/ 
CA68-10389) of EV-B81 was isolated in the USA in 1968". Subsequently, several other EV-B81 strains were 
isolated from AFP patients or healthy individuals during AFP case surveillance in China', Bangladesh'", India'", 
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Table 1 | The nucleotide sequence and deduced amino acid sequence identities between two Tibetan enterovirus B81 (EV-B81) strains 
(99279 and 99298c) and the EV-B81 prototype strain and other prototype strains belonging to enterovirus B (EV-B) 

% nucleotide identity (% amino acid identity) 
Region 

EV-B81 strain 99279 EV-B81 strain 99298c 

Prototype of EV-B8 1 Prototypes of other EV-B Prototype of EV-B8 1 Prototypes of other EV-B 



5' UTR 


78.3 


66.3-89.4 


78.4 


66.3-89.3 


VP4 


78.2 (95.6) 


67.6-79.7(71.0-82.6) 


78.2 (95.6) 


67.6-79.7 (71.0-82.6) 


VP2 


79.6 (96.5) 


64.5-73.4 (73.3-84.2) 


79.6 (96.5) 


64.3-73.5 (73.3-83.9) 


VP3 


79.8 (100.0) 


62.1-73.0 (67.5-87.4) 


79.8 (100.0) 


62.0-72.8 (67.5-87.4) 


VPl 


79.1 (94.7) 


53.8-68.8 (53.3-79.5) 


79.2 (95.4) 


53.9-69.2 (53.6-79.5) 


2A 


79.3 (93.3) 


75.1-82.8 (85.3-96.6) 


79.3 (93.3) 


75.1-82.8 (86.0-97.3) 


2B 


77.5 (94.0) 


72.6-79.8 (91.0-96.0) 


77.2 (94.0) 


72.6-79.8 (91.0-96.0) 


2C 


80.9 (96.9) 


78.0-85.7(95.1-98.1) 


81.1 (97.2) 


78.2-85.8 (95.4-98.4) 


3A 


75.2 (92.1) 


73.7-90.2 (88.7-97.7) 


75.2 (92.1) 


73.7-90.2 (88.7-97.7) 


3B 


77.2 (100.0) 


72.7-89.3 (90.9-100.0) 


77.2 (100.0) 


72.7-89.3 (90.9-100.0) 


3C 


77.9 (95.6) 


75.0-88.7(92.8-99.4) 


78.1 (95.6) 


75.2-88.7 (92.8-99.4) 


3D 


79.4 (96.7) 


77.8-86.2 (94.5-98.0) 


79.6 (96.7) 


77.9-86.1 (94.5-98.0) 


3' UTR 


84.6 


75.9-90.2 


84.6 


75.9-90.2 



Gabon'^, and Cameroon'''. To date, the only full-length genome 
sequence available for EV-B81 has been that of the prototype strain. 
Besides the prototype strain, only one entire VPl sequence and sev- 
eral partial VPl sequences of this EV type were available in the 
GenBank database. In this study, we analyzed the full-length genome 
sequences of two strains of EV-B8 1 isolated in the Tibet Autonomous 
Region of China. 

Results 

Serotyping and molecular typing of the Tibetan isolates. The 

Tibetan isolates (strain 99279/XZ/CHN/1999 and strain 99298c/ 
XZ/CHN/1999, hereafter referred to as 99279 and 99298c, respect- 
ively) were initially characterized using a standard pool of EV typing 
antisera (RIVM, the Netherlands) distributed by the World Health 
Organization. However, neither of the isolates could be neutralized 
by any of the antisera (data not shown). Therefore, the isolates were 
initially identified as "untypeable" non-polio EVs. The VPl capsid- 
coding regions of the two Tibetan isolates were then partially 
sequenced using molecular typing methods and analyzed by an 
online enterovirus genotyping tool". Both isolates were identified 
as EV-B81. 

Full-length genomic characterization of Chinese EV-B81 strains. 

The full-length genomes of the two EV-B81 strains were sequenced. 
Both were 74 1 7 nucleotides in length, encoding a polypeptide of 2 1 9 1 
amino acids. The coding sequences were flanked by a non-coding 
5' UTR of 741 nucleotides and a non-coding 3' UTR of 100 nucleotides 
followed by a poly (A) tail composed of a long sequence of adenine 
nucleotides. Alignment of the fuU-length genomes of the two Tibetan 
EV-B81 strains with the genome of the EV-B81 prototype strain 
(USA/CA68-10389) showed that they all had the same genomic 
organization and coUinear order of genomic regions. However, in 
the 5' UTR, Tibet EV-B81 strains contained three nucleotide 
insertions at positions 101, 102, and 118 and a nucleotide deletion 
at position 179. In the 3' UTR, they contained a nucleotide insertion 
at position 7328 and a nucleotide deletion at position 7341. The 
overall base compositions of strains 99279 and 99298c were 
28.03% and 28.11% A, 24.36% and 24.35% G, 23.82% and 23.73% 
C, and 23.79% and 23.81% U, respectively. The polypeptide cleavage 
sites were predicted based on the fuU-length genome sequence of the 
EV-B81 prototype strain. Table 1 shows the nucleotide sequence and 
deduced amino acid sequence identities between the Tibetan EV-B81 
strains and the EV-B81 prototype strain and other prototype strains 
within the EV-B species. The nucleotide sequences available for the 
distinct EV-B81 strains differ considerably. The complete genome 



nucleotide sequence similarity between these two Tibetan EV-B81 
strains is 99.6%, and they displayed 79.1% and 79.2% nucleotide 
identity and 94.7% and 95.2% amino acid identity with the 
prototype EV-B81 strain, respectively. 

Phylogenetic analysis of the Chinese EV-B81 strains and other 
EV-B genomes. Phylogenetic trees were generated from the 252- 
nucleotide (nucleotide 2567-2818) partial VPl coding region of 
the two Tibetan EV-B81 strains and eight other EV-B81 strains 
available in the GenBank database (Fig. 1). The two Tibetan strains 
displayed great genetic distance to another Chinese strain (strain 
142-98-YN/CHN/1998 isolated in the Yunnan province in 1998) 
and clustered with the strain BAN-10491 isolated from Bangla- 
desh. Three strains isolated from India clustered together, and 
most branches displayed long genetic distances. 

To investigate the genetic relationship between the Tibetan EV- 
B81 strains, the EV-B81 prototype strain, and other EV-B prototype 
strains available in the GenBank database, we constructed phylogen- 
etic trees based on the VPl , PI , P2, and P3 regions of the genome. The 
phylogenetic tree based on the VPl region also contained an EV-B81 
strain (N-428/IND/2008) isolated from India, for which the entire 
VPl sequence was available in the GenBank database (Fig. 2). In the 
VPl and Pi capsid regions, the two Tibetan EV-B81 strains clustered 
together with the EV-B81 prototype strain and the Indian EV-B81 
strain, confirming the preliminary molecular typing results. In the 
VPl region, the nucleotide identity between the Tibetan strains 

98j-# 99279/XZ/CHN/1999 
^IJ 99298C/XZ/CHN/1999 
I — BAN-1 049 1/BAN/AY9 19475 

I N-428/IND/2008-JN204080 

93 "-j [ NIV098601LV77/IND/2009-JX476196 

IOqI— N I V092359 1 LV1 74/1 N D/2009- JX476 1 95 

I 14622/BAN/2007-JX538042 

- I g08-024/GAB/2008-JX437661 

I 142-98/YN/CHN/1998-AB268272 

I ♦ CA68-10389/USA-AY843299 

0.02 

Figure 1 | Phylogenetic relationships based on partial VPl region 
sequences of enterovirus B81 (EV-B81). Two Tibetan EV-B81 strains 
isolated in this study (indicated by circles) and other EV-B81 strains 
(available in the GenBank database) were analyzed based on the 252- 
nucleotide (nucleotide 2567-2818) partial VPl coding region sequence. 
The strain indicated by a diamond is the EV-B8 1 prototype strain. 



SCIENTIFIC REPORTS | 4:6035 | DOI: 1 0. 1 038/srep06035 



2 



E6-DAmori 
E25-JV-4 
- E29-JV-10 




EV-B73-CA55-1988 
EV-B83-.USA-CA76-1 0392 



(b) 



'E21-Farina 
E30-Bastianni 



EV-B75-USA-OK85-1 0362 

- E1-Farouk-VP1 

— E4-Pesacek 
E20-JV-1 
E33-Toluca-3 
E13-Del Carmen 

EV-B69-Toluca-1 

CVB2-Ohio 

CVB4-JVB 
CVB1-Conn-5 
CVB3-Nancy 



(c) 



99 1 — 



100 1 

861 r 

99'— 



— E4-Pesacek-AY302557 

— E5-Nqyce-AF083069 

— E18-Metcalf-AF317694 

— E20-JV-1-AY302546 

CVB5-Faulkner-AF114383 

E19-Burke-AY302544 

— E11 -Gregory-X80059 

— E32-PR-10-AY302555 

— E24-DeCamp-AY302548 

— EV-B69-Toruca-1-AY302560 
E7-Wallace-AY302559 

— E33-Toluca-3-AY302556 

— E2-Cornelis-AY302545 

— EV-B73-CA55-1988-AF241359 

— E21-Farina-AY302547 

— E31-Caldwell-AY302554 
_ CVA9-Griggs-D00627 

— CVB3-Nancy-M16572 

— E29-JV-10-AY302552 



ri 

f)f)' — 



CVBS-Faulkner 
CVB6-Schmitt 
E24-DeCamp 
E3-Morrisey 
E12-Travis 

CVA9-Griggs 
EV-B1 07-Thra4-0349 
EV-B85-BAN00-1 0353 
E1 1-Greqory 
E19-Burke 
E7-WaMace 
E32-PR-10 
EV-B101-CIV03-10361 
EV-B111-Q0011 

99279/XZ/CHN/1999 
{• 99298C/XZ/CHN/1999 
-A N-428/IND/2008 
AE V-B8 1 - U S A/CA68- 10389 
? V-B87-B AN 0 1 - 1 0396 
EV-B88-BAN01 -1 0398 
— EV-B98-T92-1499 
E27-Bacon 
EV-B97-BAN99-1 0355 
E17-CHHE-29 
EV-B84-CIV2003-1 0603 
EV-B1 00-BAN2000-1 0500 
EV-B74- US A-CA75- 10213 
E18-Metcalf 

EV-B80-USA-CA67-1 0387 
E2-Cornelis 
E15-CH 96-51 

EV-B77-USA-TX97-1 0394 
EV-B79-USA-CA79-1 0384 
E26-Coronel-VP1 
EV-B78-W1 37-1 26-99 
EV-B86-BAN00-1 0354 
_9-Barty 
EV-B82-USA-CA64-1 0390 
E31-CaldweM 

- E14-TOW 

- E5-Noyce 

- E16-Harrington 



_ CVB1-Conn-5-M16560 

— CVB3-Nancy-M16572 

— CVB5-Faulkner-AF114383 

- CVB6-Schmitt-AF1 05342 
CVB2-Ohio-1-AF085363 

- CVB4-JVB-X05690 
. E24-DeCamp-AY302548 



E33-Toluca-3-AY302556 



. EV-B75-USA-OK85-10362-AY556070 

— E20-JV-1-AY302546 

— E1-Farouk-AF029859 

— E4-Pesacek-AY302557 

— E13-Del Carmen-AY302539 
. EV-B69-Toluca-1-AY302560 

— E21-Farina-AY302547 

— E30-Bastianni-AF162711 

— E25-JV-4-AY302549 

— E6-DAmori-AY302558 

— E29-JV-10-AY302552 



. EV-B73-CA55-1988-AF241359 
- EV-B83-USA-CA76-10392-AY843301e 



94 I — 

ft 




- CVB2-Ohio-1-AF085363 

— E15-CH 96-51 -AY302541 

— EV-B83-USA-CA76-10392-AY843301 

— EV-B79-USA-CA79-10384-AY843297 

— EV-B82-USA-CA64-10390-AY843300 

— EV-B80-USA-CA67-10387-AY843298 
— ♦ EV-B81 -USA-CA68-1 0389-AY843299 

— EV-B77-USA-TX97-1 0394-AY843302 

— E3-Morrisey-AY302553 

— E6-DAmori-AY302558 

— CVB1-Conn-5-M16560 

— E14-TOW-AY302540 
. E17-CHHE-29-AY302543 

— E16-Harrington-AY302542 

— CVB6-Schmitt-AF1 05342 

— E12-Travis-X79047 

— E25-JV-4-AY302549 

— CVB4-JVB-X05690 

— E13-Del Carmen-AY302539 

— E26-Coronel-AY302550 

— EV-B84-CIV2003-10603-DQ902712 
EV-B101-CIV03-10361-AY843308 

— E1-Farouk-AF029859 

— E27-Bacon-AY302551 
E9-Barty-X92886 

— EV-B85-BAN00-10353-AY843303 
EV-B100-BAN2000-10500-DQ902713 

— EV-B107-TN94-0349-AB426609 
EV-B87-BAN01-10396-AY843305 

— E30-Bastianni-AF1 62711 

— EV-B98-T92-1499-AB426608 

— EV-B74-USA-CA75-10213-AY556057 

— EV-B75-USA-OK85-10362-AY556070 

— EV-B88-BAN01-10398-AY843306 
EV-B97-BAN99-10355-AY843307 

— EV-B86-BAN00-10354-AY843304 

— EV-B111-Q0011/XZ/CHN/2000-KF312882 

r* 99279/XZ/CHN/1999 

mi» 99298C/XZ/CHN/1999 




. EV-B107-"fl?94-0349-AB426609 

- EV-B85-BAN00-10353-AY843303 
E3-Morrisey-AY302553 
E12-Travis-X79047 

- El 1-Gregory-X80059 

- E19-Burl<e-AY302544 

- E7-Wallace-AY302559 

- E32-PR-10-AY302555 
-EV-B101-CIV03-10361 -AY843308 

- EV-B1 1 1-Q001 1/XZ/CHN/2000-KF312882 
99279/XZ/CHN/1999 

{• 99298C/XZ/CHN/1999 

EV-B81 -USA-CA68-1 0389-AY843299 
EV-B87-BAN01-10396-AY843305 
EV-B88-BAN01 -1 0398-AY843306 
EV-B98-T92-1 499-AB426608 
E27-Bacon-AY302551 
E V-B97-BAN99- 1 0355-AY843307 
E 1 7-CH HE-29-AY302543 
EV-B84-CIV2003-1 0603-DQ90271 2 
EV-B1 00-BAN2000-1 0500-DQ90271 3 
EV-B74-USA-CA75-1 021 3-AY556057 
E18-Metcalf-AF317694 
EV-B80-USA-CA67- 1 0387-AY843298 
E2-Cornelis-AY302545 
E15-CH 96-51 -AY302541 

EV-B77-USA-TX97-1 0394-AY843302 
EV-B79-USA-CA79-1 0384-AY843297 
E26-Coronel-AY302550 
EV-B86-BAN00-1 0354-AY843304 
E9-Barty-X92886 
EV-B82-USA-CA64-1 0390-AY843300 
E14-TOW-AY302540 
El 6-Harrington-AY302542 
E5-Noyce-AT083069 

E31-Caldwell-AY302554 



E2-Cornelis-AY302545 
E 1 9-Bu rke-AY302544 
E24-DeCamp-AY302548 
E31 -CaldweirAY302554 
E32-PR-1 0-AY302555 
E17-CHHE-29-AY302543 
E29-J V- 1 0-AY302552 



E5-Noyce-AF083069 
E20-J\/l-AY302546 

EV-B69-Toluca-1 -AY302560 
. E33-Toluca-3-AY302556 
- EV-B73-CA55-1988-AF241359 

- CVB3-Nancy-M16572 
. CVB2-Ohio-1-AF085363 

- E7-WaMace-AY302559 
. El 1-Gregory-X80059 
. CVB5-Fautfcner-AF1 14383 
E15-CH 96-51 -AY302541 
CVB1-Conn-5-M16560 
E4-Pesacek-AY302557 
E18-Metcalf-AF317694 
CVA9-G rigas-D00627 
E21-Farina-AY302547 
EV-B80-USA-CA67-1 0387-AY843298 
♦_EV-B81 -USA-CA68-1 0389-AY843299 
- EV-B82-USA-CA64-1 0390-AY843300 
EV-B83-USA-CA76-1 0392-AY843301 
EV-B79-USA-CA79-1 0384-AY843297 
EV-B77-USA-TX97-1 0394-AY843302 
E14-TOW-AY302540 
E3-Morrisey-AY302553 
E6-DAmori-AY302558 
CVB4-JVB-X05690 
CVB6-Sch mitt-AF 1 05342 
E25-JV-4-AY302549 
E13-Del Carmen-AY302539 
E26-Coronel-AY302550 
E27-Bacon-AY302551 
E16-Hamngton-AY302542 
E1-Farouk-AF029859 
EV-B1 01 -CIV03-1 0361 -AY843308 
EV-B84-CIV2003-1 0603-DQ90271 2 
E9-Barty-X92886 
E12-Travis-X79047 

— E30-Bastianni-AF162711 
EV-B75-USA-OK85-1 0362-AY556070 



EV-B98-T92-1 499-AB426608 

- EV-B74-USA-CA75-1 021 3-AY556057 
99279/XZ/CHN/1999 

{• 99298C/XZ/CHN/1999 
EV-B86-BAN00-1 0354-AY843304 

- EV-B97-BAN99-10355-AY843307 
EV-B1 1 1-Q001 1/XZ/CHN/2000-KF312882 
EV-B87-BAN01-10396-AY843305 
EV-B85-BAN00-1 0353-AY843303 
EV-B1 07-TN94-0349-AB426609 
EV-B88-BAN01 -1 0398-AY843306 
EV-B1 00-BAN2000-1 0500-DQ90271 3 



Figure 2 | Phylogenetic relationships based on the VPl, PI, P2, and P3 genome regions of enterovirus B (EV-B). Two Tibetan EV-B8 1 strains (indicated 
by solid circles) and 55 other EV-B prototype strains were analyzed by nucleotide sequence alignment using the Neighbor-Joining algorithms 
implemented in the MEGA 5.0 program. Numbers at the nodes indicate bootstrap support for that node (percent of 1000 bootstrap replicates) . The open 
triangle indicates the India EV-B8 1 which has the entire VPl sequence in the GenBank database, and the solid diamond indicates EV-B8 1 prototype strain. 
The scale bars represent the genetic distance. All panels have the same scale, (a) VPl coding sequences; (b) PI coding sequences; (c) P2 coding sequences; 
and (d) P3 coding sequences. 
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99279 and 99298c and the Indian strain was 88.8% and 89.1%, 
respectively. However, in the non-capsid regions, the phylogenetic 
trees differed greatly from those in the capsid regions. In the P2 and 
P3 regions, the two Tibetan EV-B81 strains shared the highest sim- 
ilarity with the prototype strains of EV-Blll and EV-B86, respect- 
ively. This surprised us because EV-Bl 1 1 is a recently reported type 
of enterovirus'"'. The phylogenetic analysis indicated that recombi- 
nations between Chinese EV-B81 strains and other EV-B serotypes 
might have occurred. 

Recombinant structure of the Chinese EV-B81 strains. Similarity 
plots and bootscanning analyses were performed to confirm the 
recombinations between the Tibetan EV-B81 strains and other 
EV-B prototype strains. Because the two Tibetan EV-B81 strains 
share high nucleotide identity in their full-length genomes (99.6%), 
only one of the strains (99279) was used as a query sequence. It was 
compared with the EV-B81 prototype strain (USA/CA68-10389) 
and other EV-B prototype strains. In the PI coding region, strain 
99279 possessed the highest similarity with the EV-B81 prototype 
strain as expected (Fig. 3). However, in the 5' UTR, P2, P3, and 3' 
UTR regions, strain 99279 was apparently not related to the EV-B81 
prototype strain, which further confirmed the occurrence of recom- 
bination in these regions. Interestingly, a relatively high similarity 
within the 3' end of the 2C region to the 5' end of the 3C region was 
found between strain 99279 and the EV-Blll prototype strain, 
which was supported by the bootscanning analysis. And the boot- 
scanning analysis also suggested the possibility of recombination of 
small-sized fragment in 5' UTR, 3D, and 3' UTR region between 
Chinese EV-B81 strains and other EV-B serotypes such as EV-B86, 
EV-B107, and EV-B87. (Fig. 3). 

Seroprevalence of EV-B81 in Tibet. Among the 50 serum samples 
surveyed, 16 were seropositive for EV-B81 (>1:8), with a total 
positive rate of 32.0% and geometric mean titers (GMTs) of 1 : 27 
among the positive sera samples. The composition ratios for the anti- 
EV-B81 antibody titers of <1 : 8, 1 : 8-1 : 64, and >1 : 64 were 68%, 
32%, and 0, respectively. Compared with seroepidemiology studies of 
other EVs in China, the positive rate and GMTs of EV-B81 are 
apparently lower than that of other EVs such as EV-A71 and 
CVA16 in the same age group (1-5 years old)""'". 

Although both Lhasa City and Shigatse Prefecture showed low 
seroprevalence rates and low titers of anti-EV-B81 antibodies, there 
were some differences between the two areas. In Lhasa City, the 
seroprevalence rate and GMTs of the positive sera samples were 
16% and 1 : 22.6, respectively, and in Shigatse Prefecture, they were 
48% and 1 : 28.5, respectively. The positive rate in Shigatse is signifi- 
cantly higher than those in Lhasa City, while there is no significant 
difference between the GMTs of the two prefectures (seroprevalence 
rate: p = 0.0153, GMTs: p = 0.521). 

Discussion 

EV-B81 is a new type of enterovirus belonging to the species EV-B. 
The prototype of EV-B81 was isolated in the USA in 1968. 
Subsequently, several other EV-B81 strains were isolated from dif- 
ferent countries in Asia and Africa, indicating global distribution of 
the serotype. Of the eight EV-B81 strains reported in the GenBank 
database to date, six were isolated from Southeast Asia, including 
from India, Bangladesh, and provinces in southwestern China, sug- 
gesting the possible circulation of EV-B81 in this area. Detailed 
information about the prototype strain, such as information about 
its host, was not indicated in the reference publications, but almost all 
the other EV-B81 strains were isolated from patients with AFP 
except for a strain from Cameroon". Hence, there may be some 
correlation between EV-B81 and AFP. However, little research has 
been done on EV-B81 worldwide, and more data are necessary to 
unveil the biological and pathological properties of EV-B81. 



The first EV-B81 strain reported in China was isolated from the 
Yunnan province in 1998, and no other case has since been reported 
in China'^. In this study, we characterized the full-length genomes of 
two EV-B81 strains isolated in 1999 from Mangkang Prefecture of 
Tibet, located in the southeast of the Qinghai- Tibet Plateau. 
Although EV-B81 strains were reported in both Tibet and the 
Yunnan province of China, they showed great genetic diversity, the 
level of genetic diversity between two Tibetan EV-B81 strains and 
Yunnan strains was 20.4% and 20.9%, respectively, suggesting that 
EV-B81 from different lineages circulated separately in these two 
regions. For the entire VPl sequences, the nucleotide identities 
between the two Tibetan EV-B81 strains 99279 and 99298c and 
the EV-B81 prototype strain were 79.1% and 79.2%, respectively, 
further indicating that E V-B8 1 strains have high nucleotide sequence 
diversity. 

Recently, increasing numbers of studies have shown that interty- 
pic recombination is common within the species EV-A, EV-B, and 
EV-C'""^^. The two EV-B81 strains characterized in this study are no 
exception. The clustering of the two strains in the phylogenetic trees 
based on the P2 and P3 regions was inconsistent with the clustering 
in the tree based on the PI region, which suggested intertypic recom- 
bination between the two EV-B81 strains and other EV-B serotypes 
(probably EV-B86 or EV-Bl 1 1) in the P2 and P3 regions. The result 
was further corroborated by similarity plot and bootscanning ana- 
lyses. Previous studies have indicated that recombination between 
different serotypes may occur when different viruses infect and rep- 
licate in the same cell, and that recombination usually occurs among 
EV serotypes within a species'"'^"'. Hence, we can assume that the two 
Tibetan EV-B81 strains may have co-circulated with other EV-B 
serotypes, especially EV-Bl 1 1, for a period before they were isolated. 
The EV-Blll prototype was also isolated in the Tibet Autonomous 
Region in the same period. However, more data are needed to define 
the exact serotype of the donor sequence. 

Most EV infections are asymptomatic and are hence an under- 
estimated epidemic. To investigate the prevalence of EV-B81 in the 
Tibet Autonomous Region, we conducted a seroepidemiology sur- 
vey. We found that the seropositive rate and GMT were higher in 
Shigatse Prefecture than in Lhasa City (p < 0.05), indicating a geo- 
graphical difference in EV-B81 infection patterns. However, the ser- 
opositive rates and titers of anti-EV-B8 1 antibodies were overall very 
low compared with other EVs prevalent in China such as EV-A71 
and CVA16"', suggesting that the extent of transmission and the 
exposure of the population to EV-B81 were limited. However, con- 
sidering the small sample sizes (n = 50) of the survey and the limited 
number of EV-B81 strains, the possibility of cross reactivity between 
EV-B81 and antibody against other closely related EVs, in particular 
members within EV-B, could not be exclude, so further research will 
be required to draw firm conclusions. 

In conclusion, we report the full-length genome sequences of two 
EV-B81 strains isolated during AFP surveillance in the Tibet 
Autonomous Region, China. Sequence analysis revealed high genetic 
diversity in the two strains compared with the EV-B81 prototype, as 
well as intertypic recombination in the non-structural protein region 
of both strains. Although little is known about the biological and 
pathogenic properties of EV-B81 because of few research in this field 
owing to a limited number of isolates, our study provides a basis for 
further study of EV-B81. 

Methods 

Sample collection. This study did not involve human participants or human 
experimentation; the only human materials used were stool samples collected from 
AFP patients or their close contacts at the instigation of the JVIinistry of Health P. R. of 
China for public health purposes, and written informed consent for the use of their 
clinical samples was obtained from their parents of all the patients involved in this 
study. This study was approved by the second session of the Ethics Review Committee 
of the National Institute for Viral Disease Control and Prevention, Chinese Center for 
Disease Control and Prevention. 



SCIENTIFIC REPORTS | 4:6035 | DOI: 1 0. 1 038/srep06035 



4 



(a) 



5'-UTR P1/CAPSID 
I 



P2 



P3 



3'-UTR 



g| VP2 I VP3 I VP1 I 2A |2B| 2C |3A|| 3C | 



3D 



Nucleotide Position 



0 



III Ill 

1000 2000 



3000 



4000 



5000 



6000 



7000 



QUERY 

99279-XZ-CHN-1999 




1,000 



2,000 



3,000 4,000 5,000 6,000 

Position 



100 
90 



1,000 



2,000 



3,000 4,000 

Position 



5,000 



6,000 



7,000 




7,000 



- CVB1-Conn-5-M16560-genome 

- CVB2-Ohio-1-AF085363-genome 

- CVB3-Nancy-M16572-genome 
CVB4-JVB-X05690-genome 

- CVB5-Faulkner-AF114383-genome 

- CVB6-Schmitt-AF105342-genome 
■= E1-Farouh-AF029859-genome 

- E2-Cornelis-AY302545-genome 

- E3-Morrisey-AY302553-genome 

- E4-Pesacek-AY302557 -genome 

- E6-DAmori-AY302558-genome 

■ E7-Wallace-AY302559-genome 

- E9-Barty-X92836-genome 

E1 3-Del_Carmen-AY302539-genome 
Ei4-Tow-AY302540-genome 

- E17-CHHE-29-AY302543-genome 

- E18-Metcalf-AF3l7694-genome 

- E20-JV-1-AY302546-genome 

- E27-Bacon-AY302551-genome 
E29-JV-10-AY302552-genome 
E30-Bastianni-AF162711-genome 

- E33-Toluca-3-AY302556-genome 
-HEV74-USA-CA75-10213-AY556057-genome 
-HEV75-USA-OK85-10362-AY556070-genome 
-HEV80-USA-CA67-10387-AY843298-genome 
-HEV81-USA-CA68-10389-AY843299-genome 
-HEV82-USA-CA64-10390-AY843300-genome 

■ HEV83-USA-CA76-10392-AY843301 -genome 
-HEV84-CIV2003-10603-DQ902712-genome 

- HEV85-BAN00-10353-AY843303-genome 
-HEV86-BAN00-10354-AY843304-genome 
-HEV87-BAN01-10396-AY843305-genome 
-HEV88-BAN01-10398-AY843306-genome 
-HEV97-BAN99-10355-AY843307-genome 

- HEV98-T92-1499-AB426608-genome 
-HEV100-BAN2000-10500-DQ902713-genome 

HEV101-CIV03-10361-AY843308-genome 

- H EV 1 07-TN94-0349-AB426609-genome 

- HEV111-O0011-XZ-CHN-2000 



Figure 3 | Recombination analyses of complete enterovirus B (EV-B) genomes, (a) Similarity plot and (b) bootscanning analysis. A sliding window of 
200 nucleotides was used, moving in 20-nucleotide steps. The Tibetan EV-B81 strain 99279/XZ/CHN/1999 was used as a query sequence (indicated 
in the upper right corner of the image). 



The two EV-B81 strains {strain 99279 and 99298c) were isolated from stool sam- 
ples from two children living in the same village in Mangkang Prefecture of the Tibet 
Autonomous Region, China. The samples were collected in 1999, during the course of 
poliovirus surveillance activities in support of the global polio eradication initiative. 
Strain 99279 was isolated from a 3-year-old boy with AFP, and strain 99298c was 
isolated from a 5-year-old girl who was an asymptomatic contact of an AFP patient. 

For a seroprevalence study of EV-B81 antibodies, 50 healthy children < 5 years of 
age were surveyed. Fifty serum samples were collected randomly in 2010, with 
informed parental consent, by the Tibet Center for Disease Control and Prevention: 
25 samples were collected in Lhasa City and 25 samples were collected in Shigatse 
Prefecture. No children had any sign of disease at the time of sample collection. 

Viral isolation and primary identification. Stool samples from AFP patients were 
collected and processed according to standard procedures recommended by the 
World Health Organization^^. The samples were then inoculated into two cell lines, 
human rhabdomyosarcoma (RD) and a mouse cell line carrying the human 
poliovirus receptor {L20B) used to observe the development of EV-like cytopathic 
effects; the virus grew only in the RD cell line. Isolates were initially characterized by a 
micro -neutralization assay using poliovirus type-specific rabbit polyclonal antisera 
and pooled horse antisera against the most frequently isolated echoviruses and 
coxsackieviruses (National Institute for Public Health and the Environment (RIVM), 
Bilthoven, The Netherlands)^^. 

Molecular typing. For molecular typing, viral RNA was extracted from the viral 
isolates using a QIAamp Viral RNA Mini Kit (Qiagen. Germany) and stored at 



— 80'"C untU use. Reverse transcription polymerase chain reaction (RT-PCR) was 
performed to amplify the VPl coding region using PrimeScript One Step RT-PCR Kit 
Ver.2 (TaKaRa, Dalian, China) with primer pairs 008 and 013^. The PCR products 
were purified using the QIAquick PCR purification kit (Qiagen, Germany) and then 
subjected to nucleotide sequencing. Sequencing was performed in both directions 
using an ABI 3130 Genetic Analyzer (Applied Biosystems, USA), and every 
nucleotide position was sequenced at least once from each strand. The sequences were 
analyzed with the Basic Local Alignment Search Tool server at the National Center for 
Biotechnology Information and the EV serotype was determined according to a 
previously described molecular typing method^. 

Neutralizing antibody detection. Neutralizing antibodies against EV-B81 were 
detected with a neutralization test using the microtechnique and the human RD cell 
line as previously described, with some modifications^^. Serum samples were 
inactivated at 56 'C for 30 min before use, and sample dilutions of 1 : 4 to 1 : 1024 were 
assayed. Virus samples (50 |J.L) with a tissue culture infective dose (TCID50) of 100 
were mixed with the appropriate serum dilution (50 |J,L) and incubated at 36^C in a 
CO2 incubator. After incubation for 7 days, the highest dilution of serum that 
protected 50% of the cultures was recorded. A serum sample was considered positive 
if the neutralization antibody level was presented at a dilution of 1 : 8, and the GMT 
was calculated. 

Full-length genomic sequencing. Two long-distance PCR amplifications were 
performed using the Superscript III One-Step PCR-PCR system (Invitrogen). 
Reaction mixtures (50 |j,L) contained template RNA (5 ]xL), reaction buffer (25 |J,L), 
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forward [0001S48 (GGGGACAAGTTTGTACAAAAAAGCAGGCTTT)'* or E490 
(TGIGTIYTITGYRTICCITGGAT)^^] and reverse [E492 (GGRTTIGTIGWYTG- 
CCA)" or 7500A (GGGGACCACTTTGTACAAGAAAGCTGGG(T)24)'1 primers 
(1.0 ng/|iL), and RT/Platinum Taq High Fidelity enzyme (5 U). cDNA synthesis and 
pre-denaturation were carried out with one cycle of 50 "C (30 min) and 94''C (2 min). 
Amplification was carried out with 40 cycles of 94''C (15 s), 60 "C (30 s), and 68^C 
(5 min), followed by a final incubation at 68''C (5 min). The primers used for 
sequencing of the full-length genome were designed by a "primer -walking" strategy. 

Phylogenetic and bioinformatics analyses. The nucleotide and deduced amino acid 
sequences of strains 99279 and 99298c were compared with those of the prototype 
EV-B strains by pairwise alignment using the MEGA program (version 5.03)^^. 
Phylogenetic trees were constructed by the neighbor- joining method implemented in 
the MEGA program using the Kimura 2-parameter model. Regions containing 
alignment gaps were omitted from the analysis. The branch lengths of the 
dendrogram were determined from the topologies of the trees and were obtained by 
majority rule consensus among 1000 bootstrap replicates. Bootstrap values greater 
than 80% were considered statistically significant for grouping. Similarity plot and 
bootscanning analyses were performed using the SimPlot program (version 3.5.1; 
Stuart Ray, Johns Hopkins University, Baltimore, MD, USA)^^. For similarity plot 
analyses, a 200-nucleotide window was moved in 20-nucleotide steps, and 
bootscanning analyses were run with the neighbor-joining method. 

Statistical analysis. The titers of neutralization antibodies were log-transformed to 
calculate the GMTs. Chi-square test was used to compare the seroprevalence rates 
between Shigatse Prefecture and Lhasa City. Mann- Whitney U test was used to 
analyze their difference of GMTs. AU titers below 1 : 8 were assumed to be 1 : 4 for 
calculation. Differences with an error probability of P < 0.05 were regarded as 
significant. All statistical analyses were performed with IBM SPSS Statistics software 
(version 19.0). 

Nucleotide sequence accession numbers. The complete genomic sequences of the 
EV-B81 strains (99279 and 99298c) described in this study were deposited in the 
GenBank database under the accession numbers KJ755189 and KJ755190, 
respectively. 
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